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Abstract 

We calculate the spin density matrix of the vector meson produced in e'^e~ anni- 
hilation at pole. We show that the data imply a significant polarization for the 
antiquark which is created in the fragmentation process of the polarized initial quark 
and combines with the fragmenting quark to form the vector meson. The direction 
of polarization is opposite to that of the fragmenting quark and the magnitude is of 
the order of 0.5. A qualitative explanation of this result based on the LUND string 
fragmentation model is given. 







Spin effects in high energy fragmentation processes have attracted much attention [1-11] 
recently since they provide not only useful information on the spin dependence of hadronic 
interaction but also a promising tool to study the spin structure of nucleon. One of the 
important issue in this connection is the spin transfer of the fragmenting quark to the 
produced hadron which contains this quark. There are two related questions here, one is 
whether the fragmenting quark retains its polarization, the other is how one should connect 
the spin of the fragmenting quark to that of the hadron which contains the quark. 

A series of papers [1-3, 6-11], both experimental and theoretical, have been published 
on this topic. These papers all concentrate on the hyperon polarization in lepton induced 
reactions since hyperon polarization can easily be obtained in experiments by measuring the 
angular distributions of its decay products. Although the available data are still far from 
accurate and enormous enough to make a definite conclusion, the comparison of the data 
and theoretical results seems to suggest thati'0 the polarization of the fragmenting quark 
is retained in the fragmentation and that the simple SU(6) wave function can be used in 
connecting the spin of the fragmenting quark to that of the produced hadron which contains 
the quark. 

It is also interesting to note that information on polarization of vector meson can also be 
obtained from angular distributions of their decay products. Such measurements have also 
been carried out by OPAL, DELPHI and ALEPH collaboration at LEP for different vector 
mesons in e^e~ annihilation at pole [12-16]. The results show clearly that the produced 
vector meson favor the helicity zero state which implies a nonzero polarizaton of the vector 
meson in the direction perpendicular to the moving direction. There are also measurements 
on the off-diagonal elements of the spin density matrix. It is then natural to ask whether the 
data can also be understood from the same starting points as those for hyperon polarization 
and what they imply for the polarization of the antiquark produced in the fragmentation 
process of a polarized quark. 

In this note, we study these questions by calculating the spin density matrix for the 
produced vector mesons by adding the spin of the fragmenting quark (or antiquark) with 
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that of the antiquark (or quark) created in the fragmentation together. With the aid of 
an event generator JETSET, we calculate the z dependence of the polarization by taking 
all different contributions into account. Here z = 2pv/\/si where pv is the momentum of 
the vector meson, ^/s is the total e+e~ center of mass energy. We now first outline the 
calculation method in the following. 

Similar to that for hyperon polarizationS, to calculate the polarization of vector mesons 
in the inclusive process e^e" — Z° ^ V+X, we need to divide the produced vector 

mesons into the following two groups and consider them separately. (Here the subscript / 
denotes the flavor of the quark.) They are (a) those which contain the initial quark g^'s or 
the initial antiquark g^'s; (b) those which don't contain the initial quark or antiquark. The 
spin density matrix p^(-2, A;_l) for the vector meson V should be given by: 

V( 1 \ ^ {n{z,k±\aj)) V( i \ , {n{z,k±\b)) y . 
f {n{z,k±)} {n{z,k±)) 

where {n{z, k±\a, f)) and {z, k±\a, f) are the average number and spin density matrix 
of vector mesons from (a); {n{z, k±\b)) and {z, k±\b) are those from (b). {n{z,k±)) = 
X]/(^(-2, k±\a, /)) + {n{z, k±\b)) is the total number of vector mesons and k± is the transverse 
momentum of the vector meson with respect to the moving direction of initial quark or 
antiquark. Here, in contrast to the case for A hyperon production, contributions from the 
decay of heavier hadrons are very small. We just treat them in the same way as those from 
(b). 

Similar to the case for hyperon production, there are many different possibilities to 
produce the vector mesons in group (b) and the polarization in each possibility can be 
different from that in the other. Hence, it is very likely that these vector mesons as a whole 
are unpolarized. We will, just as we did in Ref.0 for hyperons, take them as unpolarized. 
This means that we simply take {z, k±\b) as a unit matrix. The spin density matrix 
p^{z, fc_|_|a, /) of the vector meson V which contains the fragmenting quark (or antiquark 
(fj) and an antiquark q (or a quark q) which is created in the fragmentation process can be 
calculated from the direct product of the spin density matrix p''/ for q'j (or p^f for (fj) and 
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that for the antiquark q (or p'^ for q). Now, we are going to formulate the calculation of 
{z^ ki_\a, /). For explicity, we will take the case V = (q^q) as an example. 
According to the Standard Model for electro-weak interaction, the initial quark or anti- 
quark produced in e~^e~ annihilation at pole is longitudinally polarized and the average 
value of polarization is Pf = —0.94 for f = d,s and b; Pf = —0.67 for f = u, Hence 
the normalized spin density matrix p'^f is given by 
\ 

(2) 



P 



C2f 



where ci/ = (l + P/)/2, C2/ = (1 — P/)/2. Here we use the helicity frame of q'j, which we 
denote by oxyz, where z-axis is taken as the moving direction of q^ in the overall center-of- 
mass frame and x-axis as the direction of the transverse momentum of the antiquark q with 
respect to the 2;-axis. The p^ is taken as the most general form in the frame oxyz, i.e.. 



( 1 



\ 



P. + lPy 1 



P. 



(3) 



where P{Px, Py, Pz) is the polarization vector of q. The spin density matrix p'^f'^ of the 
system is obtained from Eq. (2) and Eq. (3), i.e. 



py 



( 







C2fP^ 



(4) 



- Q - - 

We note that the p'^i'^ obtained in this way is in the basis of |s'^-f,s/; s"^, s|), where s'^f 
and s'^ are the spins of q^ and g, and s/ and are their z components. To obtain the spin 
density matrix for the meson, we need to transform it to the coupled basis \s,Sz), where 
s = + s^. The bases of these two representations are related to each other by a unitary 
matrix U, 

^10 0^ 



U 



71 71 
1 

yO 71 -75 0/ 



(5) 



The spin density matrix in the coupled basis is obtained from p"* = U p^f'^U , i.e. 



P' 



£1/ 
V2 



{Px ^Py 



\ 



(6) 



'^(Px+iPy) U^-PfP.) f^{Px-iPy) \{Pf-P.) 

^^(p^+iPy) C2/(l-P.) -'^{Px+iPy] 

^'^{P.+^Py) \{Pf-Pz) -'f^{P.-lPy) \{l-PfP.) J 

Hence the spin density matrix p'^{a, /) for the vector meson V in the basis of , s^) can 
be read out from Eq. (|) as follows: 



3 + PfP. 



Cl/(1+P.) '^{Px-iPy) 

'-^{Px+iPy) \{l-PfP.) '-^{Px-iPy) 

:^(p^ + iPy) C2/(l-P.) 



(7) 



In addition, we obtain also the ratio P/V of pseudoscalar meson to vector meson from Eq. 



as 



P/V 



1 - PfP. 



(8) 



3 + PfP, 

We can see that P/V = 1/3 at = and P/V > 1/3 at PjP, < 0. 

To compare p'^{a, f) with the data [12-16], we need to transform it to the helicity basis 
of the vector meson, i.e., the helicity beam frame which we denote by OXYZ. In this frame, 
Z-axis is taken as the moving direction of the vector meson and F-axis is taken as the vector 
products of Z and the e~ beam direction (see Fig. 1). This frame transformation can be 
carried out by two successive rotations: first a rotation of angle (3 around y axis which 
transforms oxyz to ox'yZ, then a rotation of angle 7 around the Z-axis which transforms 
ox'yZ to OXYZ, i.e.. 



p''{z,k^\a,f) = Dt(/3,7)p'^(a,/)D(/5,7), 

where D{l3,'~f) is the rotation matrix 
/ 



(9) 



^^cos^f sin /5 e^-^sin^l 



^^sinP cos/? —-^e^^ sin P 
e-*^ sin^ f ^ sin /? e'' 



cos^f 



(10) 



Here /3 is the angle between the momentum of the initial quark q'j and that of the vector 
meson and 7 is the angle between the y-axis and y-axis (see Fig. 1). Clearly, f3, 7, z and 
k± are related to each other by sin/3 = 2k±/{zy/s) and 
kz cos sin 6 — k± cos 6 

C0S7 = , , (11) 

' kj_ sin^ + {kz sin 6 — kj_ cos cos Oy 



where kz = y z'^s/A — k'j_ and sin 7 < for < vr and sin 7 > for > vr; 6' is the angle 
between the moving direction of the initial quark and the e~ beam direction in the overall 
center-of-mass frame and is the azimuthal angle of k± with respect to the plane of the 
initial quark and e~ beam. 

After some straightforward algebra, we obtain the following expressions for the matrix 
elements which have been measured experimentally. 

p^,{z,k^\a,f) = ^ ^ p p [1 - Pf{P. cos 2f3 + sin 2f3)], (12) 
RepX_,{z,k^\aJ) = /"f 'p p [cos 2^{Pz sin P - P.cos P) - Py sin 27] , (13) 

6 + rfTz 

Imp^_,iz,k^\aJ) = 'p ^ [sin 27(P, sin /? - P. cos /3) + Py cos 27] , (14) 

6 + rfTz 

Pe[pl^o(^) fc±|a, /) - pj„i(2:, /c±|a, /)] 

[cos 7(P^ cos 2/? - Pz sin 2/3) + Py cos /3 sin 7] , (15) 



^Pf 



3 + Pf P. 



■[sin 7(Px cos 2/3 — P^ sin 2/?) — P^ cos /3 cos 7]. (16) 



3 + P/P. 

From these results, we see already the following qualitative conclusions: 
(1) From Eqs. (8) and (|12|), we see that not only p^(2, fc_L|a, /) but also PjV depends 
on the polarization of the fragmenting quark and that of the antiquark q created in the 
fragmentation process. There is a simple relation between them, i.e.. 



p^zMaJ) = {P/V) + p""^ (P. sin/3 - P. cos/3). (17) 

We see that the relation is in general dependent on the momentum of the produced vector 
meson. If we take as usual k±= 0.35 GeV, we obtain sin/? = 2k^/ {z^) ~ 0.0077/^; for 
^/s = dlGeV'"^, which is much less than 1 for large z. Hence we have approximately 

pUz,k^\a,f) ^ P/V. (18) 



For Pz = Px = 0, we have Pqo{z, k±\a, f) = P/V = 1/3. This is the result that is expectedii 
in the unpolarized case. 

(2) From Eqs. (13-16), we see that the non-diagonal elements pY_i and {pY^ — Pq_i) 
depend not only on (3 but also on 7 which is a function of and other variables. Since for 
a given k± and a given 2;, is distributed uniformly, we should average over for these 
quantities. For k± <^ {z^/s/2), we expand cos 7 and sin 7 in terms of k±/kz and keep only 
terms up to k±/kz, then average over 0, we obtain (C0S7) ~ — tan/?cot 6'/2; (sin 7) ^ 0; 
(cos 27) ^ — tan^/3 and (sin 27) ^ 0. We insert them into Eqs. (13-16) and obtain 

RepY_,{z,k^\aJ) ^ -ii^^tan2/3(P,sin/3-P,cos/3), (19) 

3 -I- -r/Pz 

Impl_,{z,k^\aJ) ^ - /I'p p ^tan^/g, (20) 

Re[pX^,{z, fcja, /) - pI_,{z, fcja, /)] ^ -^^|^|^^^|^(P,. cos 2/? - P, sin2/?), (21) 
T \ V t 1 \ r\ y / ; I i-M y2P/Py sin /3 cot 6' 

Im[p%{z, k^\a, f) - Po,-i(^' /)] ~ 2(3 + PfP ) ' ^ ^ 

We see that they contain at least one sin/3 which is of the order of 0.01 at large z. This 
shows that these non-diagonal elements obtained in this way are all very small. More 
precisely, both Pepi _i and Jmpi _i should be smaller than 10~^ and Pe[pi — Po,-i] and 
-^"^[Pi,o ~Po,-i] should be small than 10^^. Significant nonzero results of these elements may 
imply significant interferencesEl between the fragmentation of and that of g^, which are 
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not taken into account here. Experimental resultslli^El from OPAL and DELPHI for Jmpi^_i, 
-Re[Pi,o ~ Po-i] and /m[pi q — Po,-i] seem to show no deviation from the above-mentioned 
expectations, there is however signature of deviation for Repi-i, but the statistics are still 
too low to make any definite conclusions. 

(3) From Eq. (0), we see that if = 0, i.e., the q is unpolarized in the moving direction 
of gp we obtain that 

p^oiz,k^\aJ)\p^=o = (1 - P/P.sin2/?)/3. (23) 

As we discussed above, sin 2/3 <^ 1 for large z, so that Pqq{z, k±\a, f)\p^=o ~ 1/3. This shows 
clearly that 7^ is a necessary condition for p^Q 7^ 1/3 at large z. Furthermore, neglecting 
terms proportional to sin/3 in Eq. (12), we obtain, 

p^z, k^\a, /) ^ (1 - P/P.)/(3 + PfP.). (24) 

This shows that Pqq{z, k±\a, f) > 1/3 if Pf and P^ have opposite sign. Both OPAL and 
DELPHI data [13-16] explicity show that Pqq > 1/3 for all the vector mesons except for u. 
This implies that 7^ and has the opposite sign as P/ in these cases. 

After we obtain the results for p^{z\a, f) and p^{z\b), we can calculate p^{z) if we know 
the average numbers {n{z\a, /)) and {n{z\b)). These average numbers are determined by the 
hadronization mechnism and should be independent of the polarization of the initial quarks. 
Hence, we can calculate them using a hadronization model which gives a good description of 
the unpolarized data for mutiparticle production in high energy reactions. Presently, such 
calculation can only be carried out using a Monte-Carlo event generator. We use the LUND 
string fragmentation modelEl implemented by JETSET§ to do this calculation. 

Using the event generator JETSET, we calculated {n{z\a, f )) and {n{z\b)). We use Eq. 
(p^) to calculate Pqq{z, k±\a, f) in which the only free parameter is P^. We first calculate poo 
for K*^ as a function of z since the corresponding data is available^. We found out that, to 
fit the data, the value of P^ has to be quite large. In Fig. 2 we show the results obtained 
by taking P^ = 0.48. We see that the data can reasonably well be fitted. 



Similiar calculations can certainly be made for other vector mesons such as p, 0, D*^, B* 
etc. But, for these mesons, only the average values of poo iii certain z ranges are available. 
It can easily be obtained from Eq. (1), such average values are given by, 

(Po^o) = E ^-^^-^Pooi-, f) + ^Po"o(&), (25) 

where p^(a, /) is given by Eq. ( ^41) and pj^(^) = 1/3; {n^ia, /)) and {n^ {h)) are the number 
of vector mesons from (a) and (b) in the corresponding z range respectively and (n^) the 
total number of vector meson in that z range. Using JETSET, we calculate (n^(a, /)) and 
i(ny{h)). After that we can determine the Pz which we need to fit the data on (pj^). We 
found out that, for most of the vector mesons, the resulting can be written as P^ = —o^Pf 
with a is common for most of the g^'s. In Table 1, we show the obtained results by choosing 
a = 0.51 as determined above from the data for K*^ . We see that the data can be fitted 
reasonably well except those for u and B*. The data of (pqq) can only be fitted by taking a 
as negative. The deviation in the case of B* may be attributed to the helicity flip of 6-quark 
caused by gluon radiation, which is negligible for light quarks. 

The results except that for uj can be understood qualitatively in the string fragmentation 
model. Here, it was shown thatlll the probability to creat a meson of mass m is proportional 
to e"*"^^/^, where x is the fractional energy of meson from the fragmenting quark and b is 
a positive parameter in LUND model. For a polarized quark with polarization Pf, if the 
spin of q is in the opposite direction as that of q'j, the resulting mesons can be a vector meson 
or a pseudoscalar meson with equal probabilities. The average mass is mi = (my + mp)/2, 
where my and mp are the masses for the vector meson and pseudoscalar meson respectively. 
If the spin of q is in the same direction as that of q^, only vector meson can be created and 
its mass is m2 = my. It is clear that ml < m^, thus the corresponding probability is large 
for the former case than that for the latter. This leads to a longitudinal polarization of q 
and the polarization is proportional to Pf, i.e, P^ = —a{t^Pf, 



We can see that the sign of Pz is indeed opposite to that of Pf. In case of K*^ which contains 
an initial c? or s we obtain af^^ = 0.3 for x = 0.3 by taking b = 0.58GeV^~^ as was used in 

jetsetS. 

In summary, we calculated the spin density matrix of the vector meson produced in e^e~ 
annihilation at pole from a direct sum of the spin of the polarized fragmenting quark 
and that of the antiquark created in the fragmentation process. The result for poo implies a 
significant polarization for the antiquark which is created in the fragmentation process and 
combines with the fragmenting quark to form the vector meson in the opposite direction as 
that of the fragmenting quark. The polarization can approximately be written as Pz = —<yPf 
and a ~ 0.5 for most of the mesons. 

We thank Li Shi- yuan, Xie Qu-bing and other members of the theoretical particle physics 
in Shandong for helpful discussions. This work was supported in part by the National Science 
Foundation of China (NSFC) and the Education Ministry of China. 
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TABLES 

Table 1. Spin density matrix element pqq for different vector mesons obtained from Eq.(25) using 
Pz = — 0.51P/. The data are taken from Refs.|12-lf:]. 



meson 


Poo 


data 


z range 


P 


U.oyo 




U.o < z < U.D 


P' 


0.428 


0.43 ± 0.05(DELPHI) 


z > 0.4 




0.405 


0.142 ± 0.114(OPAL) 


0.3 < z < 0.6 




0.504 


0.46 ± 0.08(DELPHI) 


z > 0.4 


</> 


0.557 


0.54 ± 0.06 ± 0.05(OPAL) 
0.55±0.10(DELPHI) 


z > 0.7 




0.415 


0.40 ± 0.02 ±0.01(OPAL) 


z > 0.5 






0.32 ± 0.04 ± 0.03(DELPHI) 




B* 


0.567 


0.33 ± 0.06 ± 0.05(ALEPH) 
0.36 ± 0.06 ± 0.07(OPAL) 


< z < 1 
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FIGURES 




Fig. 1. Illustrating graph showing the relation between the helicity frame oxyz of the initial 
quark q^^ and the helicity beam frame OXYZ of the produced vector meson V. There are three 
planes involved here: the horizontal plane determined by the incoming e^e~ and outgoing <ff<ff, 
the X — z plane, i.e., V — q^^ plane determined by the moving direction of and that of V and 
Cj/ = 6*2; X ek^/\ez X ejfcj^l, the X — Z plane determined by the moving direction of and that of 
the e~ beam and ey = x ep^_ /|ez x ep^_ \ . 
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"■3 Fmll 
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Fig. 2. Spin density matrix element poo of K*^ produced in decay obtained from Eq. (24) by 
taking Pz = 0.48 as a function of the momentum fraction z. The dotted straight Une corresponds 
to the unpolarized case (1/3); the data are from Ref. [16]. 
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Figure caption: 

Fig. 1: Illustrating graph showing the relation between the helicity frame oxyz of the initial 
quark and the helicity beam frame OXYZ of the produced vector meson V . There 
are three planes involved here: the horizontal plane determined by the incoming e^e~ and 
outgoing q^f(ff, the x — z plane, i.e., V — qj plane determined by the moving direction of 
and that of V and Cy = Cz x ekj_/\€'z x Cfcxl; the X — Z plane determined by the moving 
direction of g° and that of the e~ beam and eV = ez x ep^_/\ez x ep^_ \. 
Fig. 2: Spin density matrix element poo of K*^ produced in Z^ decay obtained from Eq. 
(24) by taking = 0.48 as a function of the momentum fraction z. The dotted straight 
line corresponds to the unpolarized case (1/3); the data are from Ref. [16]. 



Table caption: 

Tab. 1: Spin density matrix element poo for different vector mesons obtained from Eq.(25) 



using Pz = —0.5 IP/. The data are taken from Refs.|T^-[T6|. 
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